Interaction current in pp — > pp~y 
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The nucleon-nucleon bremsstrahlung reaction is investigated based on a fully gauge-invariant rel- 
ativistic meson-exchange model approach. To account consistently for the complicated part of the 
interaction current (which at present is too demanding to be calculated explicitly), a generalized 
contact current is introduced following the approach of H. Haberzettl, K. Nakayama, and S. Krewald 
[Phys. Rev. C 74, 045202 (2006)]. The contact interaction current is constructed phenomenologi- 
cally such that the resulting full bremsstrahlung amplitude satisfies the generalized Ward-Takahashi 
identity. The formalism is applied to describe the high-precision proton-proton bremsstrahlung data 
at 190 MeV obtained at KVI [H. Huisman et al, Phys. Rev. C65, 031001(R) (2002)]. The present 
results show good agreement with the data, thus removing the long-standing discrepancy between 
the theoretical predictions and experimental data. The present investigation, therefore, points to 
the importance of properly taking into account the interaction current for this reaction. 

PACS numbers: 25.10.+S, 25.40.-h, 25.20.Lj, 13.75.Cs [Phys. Rev. C80, 051001(R) (2009)] 



The nucleon-nucleon (NN) bremsstrahlung reaction 
had been studied extensively in the past mainly to learn 
about off-shell properties of the NN interaction. It 
should be clear, however, that off-shell effects are model- 
dependent and therefore are meaningless quantities for 
comparison. In fact, Fearing and Scherer [l| have shown 
explicitly for NN bremsstrahlung and related processes 
that in field theories off-shell effects cannot be measured. 

Even though the original motivation for investigating 
the NN bremsstrahlung reaction has fallen away, under- 
standing the dynamics of the NN bremsstrahlung re- 
action, nevertheless, is of extreme importance in gen- 
eral for it is one of the most fundamental processes 
involving both electromagnetic and hadronic interac- 
tions. Its importance is all the more emphasized by the 
fact that, so far, none of the existing models of NN 
bremsstrahlung can describe the high-precision proton- 
proton bremsstrahlung data from KVI [2, H[ for coplanar 
geometries involving small proton scattering angles. This 
is illustrated in Figs. QJa) an d (b), where the KVI data 
[2j for cross sections and analyzing powers are compared 
to the results of the microscopic calculations of Marti- 
nus et al. [3] and of Herrmann et al. Q. Also shown 
in Fig. Q] are the TRIUMF data for the cross section @ . 
As one can see, all the model calculations overestimate 
the measured cross sections, especially, for asymmetric 
proton scattering angles (0\ ^Q-i). In Ref. [2[, some soft- 
photon model results are also compared with the data. 
In contrast to the microscopic models, these soft-photon 
models reproduce well the measured cross-section data. 
For the analyzing powers, however, it is the microscopic 
models that describe the data much better than the soft- 
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photon models. (See also Ref. p]|, where new data from 
KVI were reported sampling a part ofphase space differ- 
ent from that for the earlier data 0, Q ; the soft-photon 
models are found to be at odds with these new data.) It 
should be noted here that, strictly speaking, the kincmat- 
ical regime of the KVI data for small proton scattering 
angles are outside of the range of applicability of Low's 
soft-photon theorem 0] . There are also a number of other 
microscopic model calculations available in the literature 
0, US ITt I JT^ L HH, 0, EH, which are dynamically similar 
to Refs. [J, y5|, addressing a variety of issues in the pp 
bremsstrahlung process. 

A detailed discussion of the status of the discrepancy 
between the theoretical and the experimental results in 
pp bremsstrahlung can be found in Refs. 0, [3 EH- This 
situation is extremely disturbing from a theoretical point 
of view, in particular, if one considers the fact that these 
data are obtained at a proton incident energy of only 190 
MeV, well below the pion-production threshold energy of 
about 280 MeV. At such a low energy, one expects the nu- 
cleonic current to be dominating by far, and baryon res- 
onances as well as meson-exchange currents should play 
only minimal roles in the reaction dynamics. 

In this work we show that the interaction current man- 
dated by gauge invariance plays a crucial role for this re- 
action and that its proper inclusion in theoretical models 
removes, to a large extent, the existing discrepancy be- 
tween the theoretical and the experimental results. For 
direct and easy comparison, our present numerical results 
— explained later in this article — are shown in Figs. [He) 
and (d), directly below those shown in Figs. [Ha) and (b) 
obtained by other authors [J, Q ■ 

The present work uses a novel approach to describe 
the NN bremsstrahlung reaction. It is derived within 
a relativistic field-theory approach by coupling the pho- 
ton everywhere possible in the underlying two-nucleon 
T-matrix determined by the corresponding NN Bethe- 
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FIG. 1: (Color online) Comparison of the KVI data @] (with corrected A y sign 0]) for pp — > ppyy in coplanar geometry at 190 
MeV proton incident energy (circles) with various theoretical predictions. The TRIUMF data [g for cross sections (triangles) 
are also displayed here. In each figure, the upper rows of panels display the cross sections and the lower rows the corresponding 
analyzing powers A y for fixed proton scattering angles, 8i and 82, and as functions of the emitted photon angle, 8 1 , in the 
laboratory frame. The two figures on the left are for the symmetric proton scattering angles, 81 =82, while the two figures on 
the right are for the asymmetric proton angles, 81 ^82- Figures (a) and (b) show the results of previous model calculations. 
The solid lines are the results of Ref. 3 including all the higher-order corrections; the dash-dotted lines pertain to the same 
model of Ref. 3 without the higher-order corrections. The dashed lines represent the results of Ref. [jj. Figures (c) and (d) 
show the same data arranged the same way compared to results obtained with the present model, i.e., the theoretical results 
shown in (c) are to be compared to those of (a) and those of (d) to (b). The dashed curves in (c) and (d) represent the nucleonic 
plus generalized contact current contributions; the dash-dotted curves correspond to the mesonic current, while the solid curves 
denote the total current contribution. Switching off the generalized four-point contact current ((3| produces the dotted curves. 
The dashed-double-dotted lines for the middle panels in (d) exhibit the parameter sensitivity explained in the text. 



Salpeter equation. The basic idea of this formalism is 
the same as what had been introduced by Haberzettl, 
Nakayama, and Krewald for pion photoproduction, 
based on the field-theoretical approach of Haberzettl [13] • 
The full bremsstrahlung amplitude can be written as 

M M = (TG + 1)J M (1 + G Q T) , (1) 

where the NN T-matrices on the left and right mediate 
the final-state (FSI) interaction and the initial-state in- 
teraction (ISI), respectively; Go denotes the intermediate 
propagation of two free nucleons. The current 

= d^GoV + VG d" + - VG d>*GoV (2) 



is the basic photon production current off the two nu- 
cleons composed, in general, of nucleonic, mesonic, and 
baryon-resonance currents, in addition to contact-type 
interaction currents, as shown in Fig. HJa) for the ex- 
ample of single-meson exchanges. The two disconnected 
nucleonic contributions subsumed in d M are shown in 
Fig. [2jb); V is the NN interaction and describes 
the photon coupling to the internal mechanisms of the 
interaction V. We emphasize that the generic structures 
of Eqs. ||T]| and ||2J| are exact for any type of NN interac- 
tion. The proof follows straightforwardly from applying 
the "gauge derivative" procedure given in the Appendix 
of Ref. to the connected part GqTGq of the NN 
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FIG. 2: (a) Basic production amplitude of eq. (J2j) for NN — > NN*y used in the present work. (Time proceeds from right 
to left.) The terms below the diagrams correspond to the respective ones in Eq. ([2]), with (b) showing the photon coupling to 
both intermediate nucleons subsumed in eP. The diagrams for the lower nucleon line analogous to diagrams 1, 2, and 4 are 
suppressed. N denotes the intermediate nucleon and M incorporates all exchanges of mesons n, n, p, ui, a, and a (former 8). 
In general, contributions to the N N interaction more complex than single-meson exchange may be considered as well. External 
legs are labeled by the four-momenta of the respective particles; the hadronic vertices s, u, and t (with the labels alluding to 
the corresponding Mandelstam variables) correspond to the same kinematic situations, respectively. The first two diagrams on 
the right-hand side describe the so-called nucleonic current and the meson-exchange current is depicted by the third diagram. 
The fourth diagram contains the NM — > Nj four-point contact current of Eq. @, labeled "c" in the diagram. The first four 
diagrams correspond to the complete gauge-invariant description for the process NM — > Ny for the upper nucleon line. The 
last diagram (labeled "C") stands for the five-point contact-type current in general necessary to preserve gauge invariance of 
the entire amplitude, including the AW ISI and FSI contributions. 



Green's function. For NN interactions based on single- 
meson exchanges, the structure of J p depicted in Fig.[2]is 
complete. We draw particular attention to the contact- 
type currents appearing in the fourth and fifth term in 
Fig. [3Ja) . In general, they possess very complex internal 
dynamical structures that, at present, cannot be taken 
into account explicitly. However, they are constrained by 
gauge invariance. The corresponding constraint for the 
four-point current in the fourth diagram is based on the 
generalized Ward-Takahashi identity for the subprocess 
NM — * Nj. 1 By contrast, the five-point contact current 
depicted in the last diagram of Fig. Ufa) is constrained by 
demanding gauge invariance of the full bremsstrahlung 
process using the fact that the subprocesses already sat- 
isfy their respective gauge-invariance constraints. 

For the current J M considered here shown in Fig.[2ja), 
in addition to the first two diagrams on the right-hand 
side with intermediate nucleons marked N, one may also 
consider contributions from intermediate baryon reso- 
nances. However, for the present application to the KVI 
data Q at 190 MeV incident proton energy, we expect 
their contributions to be minimal, and we therefore have 
omitted such contributions. 2 



1 See the corresponding discussion in Ref. [id for the case of 
pion photoproduction that exemplifies the structure of the 
four-point contact current appearing in the fourth diagram of 
Fig. [2j a )- This case is relevant here because the dynamics of 
NN bremsstrahlung can be largely understood as a meson cap- 
ture process where the captured meson originates from a specta- 
tor nucleon, as can be seen in the first four diagrams in Fig,l2Ta), 

2 For the A, in particular, we point to the results of Ref. 4] shown 
in Figs. [TJa) and (b) with (solid curves) and without (dashed- 
dotted curves) higher-order effects, where the former includes 
the A as well as the mesonic current contributions. They show 
that the A would have only a minor effect in the low-energy 
regime investigated here. This therefore cannot affect our overall 
conclusions. At higher energies, effects of the A resonance have 
been investigated by de Jong et al. [lC| . 



Our present approach is fully relativistic, employing 
a covariant three-dimensional reduction of the Bethe- 
Salpeter equation underlying the NN interaction 
As an immediate consequence of this reduction, one 
finds that the five-point contact current [last diagram in 
Fig. Ufa)] must not contribute to the four-divergence of 
the full amplitude, i.e., it must be fully transverse. (This 
is not true for the amplitude M M evaluated in a full four- 
dimensional framework. The five-point contact current 
then is essential to maintaining gauge invariance.) In 
our present calculation, therefore, because it has no bear- 
ing on gauge invariance, we have dropped the five-point 
contact current entirely. Further details of the present 
formalism will be reported elsewhere. 

In view of this finding, gauge-invariance constraints 
come to bear only on the four-point interaction current 
describing the NM — > N*y subprocess in the fourth di- 
agram of Fig. HJa). To avoid having to deal with its 
very complex microscopic dynamical structures [l6j], we 
employ a generalized contact term that is constructed 
such that the resulting full amplitude satisfies the gener- 
alized Ward-Takahashi identity necessary to ensures full 
gauge invariance. The details of the contact term em- 
ployed here are discussed later in this article in connec- 
tion with Eq. ([31). The only mesonic currents [cf. third 
diagram in Fig. H[a)] contributing to pp bremsstrahlung 
are those arising from the anomalous couplings which 
cannot be obtained from coupling the photon to the un- 
derlying NN interaction. Following Ref. [l6[, we include 
the 7rp7- and 7rw7-exchange contributions in the present 
work which are the dominant mesonic currents for pp 
bremsstrahlung. 

In the present work, we use the OBEP-B version of 
the Bonn AW interaction with its parameters slightly 
readjusted to reproduce the low-energy pp scattering 
length (see also Ref. (H). This interaction contains 
only nucleon and meson degrees of freedom. The rea- 
son for choosing this interaction is that, apart from the 
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phenomenological form factors at the nucleon-nucleon- 
meson (NNM) vertices, it is a fully microscopic meson- 
exchange model and, as such, the photon can be attached 
consistently and uniquely to every part of the NN inter- 
action except, of course, to the form factors. The latter 
mechanism introduces necessarily an ambiguity in how 
the photon couples to this interacting system. As alluded 
to previously, we account for it through the generalized 
four-point contact current given below. We note once 
more that with such a simple explicit microscopic model 
of the NN interaction based on a three-dimensional re- 
duction of the underlying Bethe-Salpeter equation, no 
five-point contact current is required to maintain gauge 
invariance, i.e., all currents can be calculated explicitly, 
with the exception of the four-point contact current aris- 
ing from the photon coupling to the NNM vertices, 
as pointed out. A five-point contact current may be 
necessary when one uses a more sophisticated NN in- 
teraction that includes two-nucleon irreducible contribu- 
tions, as would be the case, for example, with interac- 
tions including explicit A degrees. For such interactions, 
and also for those based on purely phenomenological ap- 
proaches where the underlying microscopic structures are 
not known, the introduction of a five-point contact cur- 
rent becomes an unavoidable procedure. 

As alluded to in footnote [TJ it is instructive to con- 
sider the present approach from a different point of view. 
Ignoring for the moment the complication of a possible 
five-point contact current, the basic photon production 
amplitude J M can be viewed as being obtained from the 
inverse of the photoproduction reaction, i.e., a meson 
capture reaction, where the captured meson originates 
from a second nucleon that is a spectator to the capture 
reaction. In this picture, the nuclcon-mcson FSI is ac- 
counted for effectively by the four-point contact current. 
This is discussed in detail in Ref. flq . In this sense, one 
may think of the full bremsstrahlung amplitude given by 
Eq. |T]) as being constructed from the two basic build- 
ing blocks, the respective amplitudes for the NN — > NN 
and the NM — > Nj reactions, analogous to the approach 
described in Ref. [2(J for calculating the NN -> NNrj 
process. Of course, in this simplified picture one ignores 
complications arising from the identity of the nucleons. 

For the details of the derivation of the generalized four- 
point contact current involving a pseudoscalar meson ex- 
change, we refer to Ref. [l6j |. Here, the required exten- 
sion — to account for the virtual nature of the incoming 
and outgoing nucleons and the exchanged meson at the 
four-point vertex — is accomplished following the work 
of Ref. [l8| . Furthermore, in the present work, the gen- 
eralized four-point contact current is extended to include 
also the scalar and vector meson exchanges. Schemat- 
ically, suppressing all Lorentz indices, we have for the 
fourth diagram in Fig. d^a) 

J £ =£(/»r 2 )iAa, [ejtf/xtr^ + riCf] , (3) 

M 

where the summation runs over the exchanged mesons 



between the interacting nucleons 1 and 2. Nucleon 2 
(described by its phenomenological form factor fa t and 
its Lorentz operator structure T 2 ) is here the spectator 
nucleon that supplies the meson (described by the propa- 
gator Am) participating in the NM — > A^7 subprocess at 
the other vertex. The corresponding four-point interac- 
tion current is described by the expression in the square 
brackets here. The first term contains the phenomeno- 
logical form factor f\t and current operator r£ , i.e., this 
describes the usual Kroll-Rudcrman-type NNM^j con- 
tact vertex. The second term contains the Lorentz oper- 
ator Ti for the NNM vertex of nucleon 1 and the contact 
current [TEG]} 



ei 



(p' x + k) 2 -p' x 

(2pi - ky 
(pi - k) 2 - p\ 

(2q - ky 
(q — k) 2 — q 2 



(fin -Ft) 



(fit ~ Ft) 



(4) 



that is necessary for the preservation of gauge invariance 
because the nucleons have structure described by form 
factors. The factors e' l5 ei, and em stand for the com- 
bined charge-isospin operators of the nucleons and me- 
son at the four-point vertex. The subtractions of F\ in 
Eq. (HJ) are necessary to render Cf pole-free; F\ is a phe- 
nomenological function chosen here as 



F 1 =R 1 - h 



(Ri ~ S s fi s ){Ri - S u f lu ){R! - S t fit) 

Ri 



, (5) 



which is manifestly crossing symmetric. For nonzero 
charges e x , one has 5 X = 1 and zero otherwise; f\ x de- 
note the hadronic form factor for the specified kinematics 
x = s, u, t [cf. Fig. [DJa)]; p[ and p\ are the four-momenta 
at the four-point vertex of the outgoing and incoming 
nucleon, respectively; and q and k stand for the four- 
momenta of the exchanged meson M and the emitted 
photon, respectively. In Eq. l[5]). 



Ri = l + e- 2/a (/-l) 



(6) 



with z = [(1 - 5 a fu)(l - <S„/iu)(l - 5 t fit)} 2 , and / = 
F(p[ 2 ,p 2 ,q 2 ) denotes the hadronic form factor with the 
momentum arguments as indicated. Of course, an ex- 
pression similar to Eq. ([3]) exists for the photon emerging 
from the NNM vertex of nucleon 2. [The latter diagram 
is not shown in Fig. [H(a).] 

The only free parameters of our model are the param- 
eters h and a appearing in Eqs. {5} and ([6]) in the gen- 
eralized four-point contact currents involving the scalar, 
pseudoscalar, and vector meson exchanges [cf. Fig. 0(a), 
fourth diagram]. In principle, these parameters may both 
be chosen independently for different exchanged mesons, 
in addition to being functions of momenta at the four- 
point vertex. In the present work, we take them to be 
constant and equal for all the exchanged mesons for the 
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sake of simplicity. Their values of h = 2.5 and a = 1000 
have been adjusted to reproduce the KVI cross section 
data 0. All coupling constants and form factors at the 
hadronic vertices are consistent with the Bonn NN inter- 
action [l9| we use here for the FSI and the ISI. The only 
exception is the NNu> coupling constant, gNNu, and the 
cutoff parameter, A ff , in the NNir form factor entering in 
the mesonic current. Following the discussion in Ref. [2l| . 
we take qnnuj = 10 and A^ = 900 MeV. Also, following 
the work of Refs. [2(| Hl| , we employ a form factor for the 
off-shell nucleon in the basic photon production current, 
J* 4 , with a cutoff parameter value of Ajy = 1000 MeV. 

In Figs.QJc) andfljd), we show the present results for 
the cross sections and analyzing power for the symmetric 
(left panels) and asymmetric (right panels) proton scat- 
tering angles. As one can see, we now reproduce well 
the cross-section data for both symmetric and asymmet- 
ric proton scattering angles (solid curves). The dotted 
curves correspond to the results when the generalized 
contact current is switched off. This illustrates the im- 
portance of taking into account the interaction current 
properly, a feature that has been ignored in all the earlier 
models. In fact, our results without the contact current 
are very similar to the ones obtained in earlier models, 
especially for asymmetric proton scattering angles [com- 
pare, in particular, the top-row panels in Figs, [lja) and 
UJb) with those in Figs.QJc) andQJd), respectively]. The 
contact current, however, does not affect the analyzing 
power significantly and thus leaves room for further im- 
provements. 

To provide some idea about how sensitive the present 
results are to the fit parameters h and a, the dashed- 
double-dotted curves in the middle panels of Fig. DJd), 
for 6>i = 8°, 2 = 19°, correspond to h = 3.0 in Eq. ©. 
Note that the agreement with the data improves for both 



the cross section and the analyzing power. The calcu- 
lated results are rather insensitive to the parameter a 
in Eq. ([5]), once it is in the correct order-of-magnitudc 
range. 

The Bonn JVJV interaction [l9( used here does not in- 
corporate Coulomb effects. To account for them, we re- 
peated the calculation using instead the Paris NN inter- 
action [22| , which includes the Coulomb interaction fully 
as described in Ref. [f| • We obtain results practically the 
same as the ones shown in Figs. [Tic) andQJd) if we read- 
just the parameter value h in Eq. ([5]) to h = 2.0. The 
only noticeable Coulomb effect is a downward bending of 
the cross section near 7 = 180° in the 0i = 02 = 8° ge- 
ometry, which is in agreement with earlier findings [5j, [9( . 

In summary, our present results essentially resolve the 
long-standing discrepancy between theoretical and ex- 
perimental results in the pp bremsstrahlung reaction. 
The new feature of the present model responsible for 
bringing the theoretical results in line with the mea- 
sured high-precision cross-section data from KVI 0] is 
a generalized four-point contact current that accounts 
for the interaction current in the NM — > V7 subpro- 
cess and that is constructed in such a way that the re- 
sulting full bremsstrahlung amplitude obeys the Ward- 
Takahashi identity and thus is gauge invariant. 

In view of its phenomenological nature we do not ex- 
pect the particular contact current employed here to 
provide a definitive resolution of all problems in NN 
bremsstrahlung processes. Nevertheless, our present re- 
sults show that the interaction current is a necessary in- 
gredient for bremsstrahlung calculations that cannot be 
neglected. 

This work is supported in part by the FFE-COSY 
Grant No. 41788390. 
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